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Abstract 



The recent analyses of the world neutrino data, including the T2K and MINOS 

1-^ ' results, point toward a statistically significant deviation of ^13 from zero. In this paper 

we present a SUSY SU(5) model based on the discrete 5*4 group which predicts a large 

Qh| 6*13 ~ 0{\c), A(7 being the Cabibbo angle. The other mixing angles in the neutrino 

(— I ■ sector are all compatible with current experimental data. In the quark sector, the 

entries of the CKM mixing matrix as well as the mass hierarchies in both up and down 

pi I quark sectors are well reproduced and only a small enhancement is needed to reproduce 

>■ Ac. 
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1 Introduction and description of the model 

The T2K experiment has recently observed six events which, after all selection criteria at 
the far detector, are a strong indication of u^ — )• Uf, flavour transition [!]. In a three flavour 
scenario with lAmjsl = 2.4 x 10^^ eV^ and maximal atmospheric mixing, these data are 
consistent with a non-vanishing 6^13 at 2.5cr, with 0.03(0.04) < sin^ 26'i3 < 0.28(0.34) for 
normal (inverted) hierarchy and best fit value at 0.11 (0.14). This result goes in the same 
direction of other, statistically less significant, analyses where a not-so-small reactor angle 
emerged from global fits to the available neutrino data [21 |3l H]. 

Hints for ^13 > have been confirmed in [5], where an analysis of global neutrino data, 
including the latest T2K and MINOS results [6], provides 

, r 0.021 ±0.007, old reactor fluxes , , 

sin^^i3=<^ (la), 1 

1 0.025 ± 0.007 , new reactor fluxes ^ ^ ^ ^ 

which corresponds to a > 3cr evidence for a non-vanishing reactor angle, for both old and 
new reactor neutrino fluxes [7]. A large value of ^13 seems to disfavor the picture where the 
tri-bimaximal (TBM) mixing pattern [8],(tan^^23 = l,tan^0i2 = | and sin0i3 = 0), is a 
good first order description of the data. In fact, models based on A4, 5*4, T7, T' and so on 
(see [9j for a review) have the common feature that, at the next level of approximation, all 
the three mixing angles receive corrections of the same order of magnitude, which is basically 
fixed by the experimentally allowed departures of 612 from its TBM value, at the level of 
0(A^), with Ac being the Cabibbo angle. As a consequence, 6^13 (and also the deviation 
of 623 from the maximal value) is expected to also be at most of (9(A^), which is now only 
marginally allowed (but see [10] for some examples of large 613 in the context of TBM). As 
it has been pointed out in [TT], a large value of the reactor angle can be achieved if the 
Bimaximal (BM) mixing option, tan^ ^23 = l,tan^^i2 = 1, sin ^13 = 0, is the correct first 
order approximation to describe the neutrino mixings; corrections from the charged lepton 
sector must be large enough, of 0{\c), to move 613 toward large values and to reconcile the 
value of 612 with the experiments but much smaller in the atmospheric sector in order not 
to destroy the agreement with maximal mixing. This pattern of subleading corrections is in 
fact realized in [TI], giving: 

sm'9i2 ~ l + 0{Xc) sm'923-^^ + 0{Xl) sin^i3 ~ 0(Ac) • (2) 

In this paper we consider a GUT extension of such a model, based on the SU(5) gauge group, 
with the aim of extending the 6*4 symmetry to describe the quark sector and maintaining at 
the same time the relations in eq.([2]). Earlier attempts to describe quarks and leptons in the 
context of SU(5)x5'4 can be found in [13] where, however, the main goal was to reproduce 
the TBM. Similar tentatives in SO(IO) and Pati-Salam have been discussed in [12] and 
[13], respectively. The model proposed here is the first attempt to get the BM pattern in 
a SU(5) context based on 6*4. In order to formulate a reahstic GUT model, we work in a 
supersimmetric scenario in 4+1 dimensions where problems related to the breaking of a grand 
unified symmetry (like the doublet-triplet splitting problem and the proton decay) can be 
efficiently solved [TS]. In the simplest setting, the fifth dimension is compactified on a circle 
S^ of radius R in such a way that the gauge fields, living in the whole 5D space-time, are 
assumed to be periodic along the extra dimension only up to a discrete parity transformation 
fl under which the gauge fields of the SU(3)xSU(2)xU(l) subgroup are periodic and possess 
a zero mode while those of the coset SU(5)/SU(3)xSU(2)xU(l) are antiperiodic and form 
a Kaluza-Klein tower starting at the mass level 1/R. For a 4D observer, these boundary 



conditions break SU(5) down to the Standard Model (SM) gauge group at a GUT scale of 
order 1/R. The doublet-triplet splitting problem is solved if the parity Vt is extended to 
the Higgs multiplets H^ and i/^-, also assumed to live in the bulk, in such a way that the 
electroweak doublets are periodic (and then have zero modes), whereas the colour triplets 
are antiperiodic (getting masses of order 1/-R). To reduce N = 2 SUSY (induced by the 
original A^ = 1 SUSY in five dimensions) down to A^ = 1 it is necessary to compactify the 
fifth dimension on the orbifold S^/Z2 rather than on the circle S^. The orbifold projection 
eliminates all the zero modes of the extra states belonging to A^ = 2 SUSY and also those 
of the fifth component of the gauge vector bosons. The zero modes we are left with are the 
4D gauge bosons of the SM, two electroweak doublets and their A^ = 1 SUSY partners. For 
the gauge vector bosons and the Higgses we adopt this setup, which is described in detail in 
refs. [TB]. For the remaining fields we have much more freedom [TCI [IZ]: they can be located 
in the bulk, at the SU(5) preserving brane y = 0, or at the SU(5) breaking brane y = ttR. 
In our construction based on 5*4, the three 5 are grouped into the 5*4 triplet F, while the 
tenplets Ti, T2 and T3 are assigned to the singlet 1 of S'4. We choose to put the tenplet of the 
first two families Ti,T2 in the bulk and all remaining A^ = 1 supermultiplets on the SU(5) 
preserving brane at y = 0. To obtain the correct zero mode spectrum with intrinsic parities 
compatible with symmetry and orbifolding, one must introduce two copies of each multiplet 
with opposite parity Q in the bulk. Therefore Ti^2 is a short notation for the copies Ti^2 and 
T[2', the zero modes of T12 are the SU(2) quark doublets Qi,2, while those of T^' 2 are 1/^2 
and E^2- This setup has a two- fold advantage: an automatic suppression of the Yukawa 
couplings for the fields living in the bulk and the breaking of the mass relation rrie = mj for 
the first two fermion familie^j. The latter is a consequence of using the copies Ti^2 and T{ 2 
for the down quarks and charged leptons, respectively, whereas the former originates from 
the fact that a bulk field B and its zero mode B^ are related by: 



B = ^=B^ + ... (3) 




where dots stand for the higher modes. This expansion produces a (geometrical) suppression 
factor 

s = -jL= < 1 (4) 

entering the Yukawa couplings depending on the field 5°. This applies also to the Higgs 
vevs; since all the matter fields (but T\^2) are localized at y = 0, what matter for the Yukawa 
couplings are the values of the vevs at y = 0: 



(^5(0)) = ^ (%(0)) = ^%. (5) 



vl ,„ .„„ V)j^ 



A similar setup (to which we refer for further details) has been used in jTH] in the context of 
a SUSY SU(5)x744. To break the S'4 symmetry, we consider a set of SU(5)-invariant flavon 
supermultiplets: three triplets ^i,^v (3i), Xt (82) and one singlet ^y. The alignment of their 
vacuum expectation values (vevs) along appropriate directions in flavour space will be the 
source of BM lepton mixing. For this to work, we employ a cyclic Z3 symmetry which allows 
the fields i^y and ^y to be the only ones responsible for neutrino masses at leading order. 



^It is interesting to observe that, to break the SU(5) relation rrie = fn^ in the case Ti,2 are also localized 
on the brane, one could introduce a 45 representation for the Higgs fields, also propagating in the bulk. 
However, it turns out that this representation not only contains zero modes for the electroweak doublet but 
also for colored Higgses, thus reintroducing a sort of doublet-triplet splitting problem. 



The GUT Higgs fields H^ and H-^ are singlets under the family symmetry but charged in 
the same way under Z3, so that they are distinguished only by their SU(5) transformation 
properties. To achieve a realistic mass spectrum, beside the geometrical suppression factors, 
we also exploit the Froggatt-Nielsen mechanism. The tenplets Ti and T2 are charged under a 
U{1)fn flavour group, spontaneously broken by the vevs of two fields and 0' both carrying 
U{1)fn charges —1 and transforming as a singlet of 5*4. The assignment of the fields under 
SU(5) and the discrete group 5*4 x Z3 is summarized in TabiU Before closing this section, it 
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Table 1: Matter assignment of the model. The symbolhr(hu) indicates that the corresponding 
fields live on the brane (bulk). 



is interesting to outline that, as long as the U(l)i(; symmetry remains unbroken, dangerous 
operators that could spoil the solution of the doublet-triplet splitting problem or mediate 
proton decay are forbidden to all orders. In fact, given the U(1)r assignments in TabU], 
the mass term H^H^ has U(l)ij = and cannot be included in the superpotential of the 
effective A^ = 1 SUSY, which should have U(l)ij charge +2, to compensate the i?-charge 
—2 coming from the Grassmann integration measure d^6. Also, all renormalizable baryon 
and lepton number violating operators, such as FH5 and FFT, are not allowed, and the 
dimension five operator FTTT, leading to proton decay, has -R-charge +4 and therefore is 
absent. The paper is organized as follows: in SectEJwe discuss the vacuum alignments of 
the fiavon fields at LO and NLO, since they are the necessary ingredients to build the mass 
matrices for charged and neutral fermions; in Sectl3]we show how to get a realistic pattern 
of the mass ratios in the down quark and charged leptons, also computing the left-handed 
rotations in both sectors, needed to build the CKM matrix and correct the neutrino mixing 
matrix. SectS] is devoted to the up- type quarks whereas in sectJS] we discuss the neutrino 
sector and some phenomenological implications of the LO and NLO mass matrices. In SectE] 
we draw our conclusions. 



2 Vacuum alignment 

We solve the vacuum alignment problem using the method first introduced in [19]. Within 
this approach, a continuous t/(l)i? symmetry is introduced, under which matter fields have 
R = +1, while Higgses and fiavon fields have R = 0. Such a symmetry will be eventually 
broken down to the R-parity by small SUSY breaking effects which can be neglected in the 
first approximation. The required vacuum alignment is obtained introducing the so-called 
driving fields with R = +2, which enter linearly into the superpotential. We use here the 
same fiavon content and S4 property transformations as in [11], that is two triplets cp'l and 
X°, one doublet ijj'^ and one singlet ^°. We assume that the family symmetry is broken at an 
energy scale where SUSY is still an exact symmetry; this allows to deduce the alignment of 
the fiavon fields from equations arising from setting the F-terms of the driving fields to zero. 



All the multiplets but the flavon ones have vanishing vevs and set to zero for the present 
discussion. We regard the U{1)fn Froggatt-Nielsen flavour symmetry as a local symmetry, 
assuming that other vector-like multiplets (not specified here) are introduced to remove the 
anomaly associated with U{1)fn- Within these assumptions the relevant part of the scalar 
potential of the model is given by the sum of the F-terms and of a D-term: 



with 



V = Vp + Vn , (6) 

dw 



Vf = E 



dipi 



2 

(7) 



2.1 Leading order 



The LO superpotential responsible for the vacuum aligment is equal to the one quoted in 

m- 

+ C [Ml + M'^i, + g:,{ipuipu) + g^i.iu] + (8) 

+ /4 {x%^ad^2) ■ 

We parametrize the triplet fiavon vevs as 

(0) =V4>\ 



In the SUSY limit, the vacuum configuration is determined by the vanishing of the derivative 
of Wd with respect to each component of the driving fields. The set of such equations for 
the minimum of the potential can be divided into two decoupled parts: one for the neutrino 
sector (involving ify and ^y) and one for the charged lepton sector (driven by iff, and x^)- 
We do not report here such equations (since they are equal to [Hj) and only quote the final 
results: 





(9) 




(e.) = v^ 



where the various f^'s obey: 
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(10) 



with v^ undetermined. The D-terni is given by: 

Vd = \{Ml, - gFN\e\'' - gFN\e'\^ + ...)' (11) 

where gp^ is the gauge couphng constant of U{l)p]^ and Mpj denotes the contribution of 
the Fayet-Ihopoulos term. We have omitted the SU(5) contribution to the D-term, whose 
vev is zero. There are SUSY minima such that Vp = Vd = 0. The vanishing of Vd requires 

gFN\e\^ + gFN\e'\'' = mIj . (12) 

If the parameter Mp^ is positive, the above condition determines a non-vanishing vev for a 
combination of 9 and 6' . 

2.2 Next to leading order 

The next level of approximation is different from [H] . The corrections to the flavon super- 
potential can be expressed as follows: 

i=l i=l i=l 1=1 

where the following operators contribute to the quartic invariants < /^ , /,■ , I^' , I^*^ > (we 
do not specify here the different contractions among the fields, see appendix A for details): 

^f" : fl, V^i xh ¥>%, <^hi^ x'^e'^e, ¥>uC 

if : ^l, e^, ifl iflU x'm (14) 

h ' ■ ^]^u, ^U^, Xl^u, X%, 'fu^lXt, iv^tXt 







The relevant feature of such corrections is the presence of terms where flavons of the neutrino 
and charged lepton sectors mix to each other. To get the NLO vacua, we parametrize the 
vev shifts according to: 

{Vi) = v^^\ 1 + 5^^ I {xi) = v^ I (5^2 I (15) 

5. 



5^3 / V 1 



'Ul 



M = v^^\ 1 + 6,, I {^^) = y^(l + Sv^) (16) 

-1 + 5,3 

where all 6i are smaller than 1 , and look for solutions by imposing the vanishing of the first 
derivative of Wd + Aw^ with respect to the driving fields. After some algebra, we get the 
following results: 
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where all the non- vanishing 5i are proportional to a combination of fj/A; this ratio also fixes 
the relative magnitude between the leading order components of the flavon vevs and the 5iS. 
To make this more transparent, we rescale the perturbations 5i — )■ e' 5i, with now 5i ~ 0{1). 
It is important to observe that the corrections to the second and third component of ip^ 
are opposite to each other and can be reabsorbed into the leading order results. The same 
happens for the non- vanishing corrections to the second component of ifi. According to what 
discussed, the fiavon vev structure that will be used to determine the fermion masses can be 
summarized as follows: 
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(20) 



-1 



The magnitude of the fiavon vevs as well as of the e' perturbations by will be discussed in 
the next section. 

3 Down quarks and charged lepton mass matrices 

The superpotential built with operators with two-fiavon insertions (beside the fiavons carry- 
ing Froggatt-Nielsen charges) allows to determine the relevant features of the down quarks 
and charged lepton mass matrices. According to the discussion in the Introduction, we use 
the same notation for the couplings involving the tentplet Ti 2, although it should be under- 
stood that they are different for down quarks and charged leptons. The superpotential in 
this sector reads as follows: 
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where we took into account that the fields living in the bulk have mass dimension 3/2. From 
the superpotential in eq. ([8]) we expect a common order of magnitude for the vev's (scaled 
by the cutoff A): 



v^^ ~ t;^ ~ e 



1! ^^ 



f ^ ~ e . 



(22) 



although, due to the different minimization conditions that determine {v^^,v^) and {v^^,v^), 
we may tolerate a moderate hierarchy between e and e. Similarly the order of magnitude of 



{9) and {9') is in principle unrelated to those of e and e. In the following, we assume that 
e = e and use the short-hand notations for the ratio among flavon vevs and A: 

^ = ^ = ^ = ? = ^ ^ = ^ ^ = ^'- (23) 

A A A A A A ^ ' 

Disregarding for the moment all 0{1) coefficients and taking the vev {H-^) = Vd, the mass 
matrix (in the ipi^^di^R convention) obtained from this Lagrangian is: 

/ se{t^ + tH' + t'3) s tt'^ s tt'{t + t') \ 
md = VdSe I set' st st , (24) 

\ e 1 J 

where s is the suppression volume factor. If e, t and t' are relatively close in magnitude, we 
can easily recognize that the down quark and charged lepton mass hierarchies are given by: 

nih : rUs : nid ~ rrir : m^ : nie ^ 1 : st : se (t^ + t^t' + 1'^) (25) 

whereas the mixing angles in the down sector can be estimated to be: 

9i^ ~ t'^ 9i^ ~ s tt'it + 1') 9i^ ~ st . (26) 

Notice that, by transposition, we also get an estimate of the charged lepton mixings, given 
by: 

9[, ~ e (^0 9[, r^e ^ ~ 0. (27) 

Then, a realistic pattern of fermion masses and mixings can be achieved requiring that 

t ~ t' ~ A , (28) 



nu }- r-u 



where A = Ac- Given that e' 5i ~ fj/A, we can also fix e' ~ A. Before discussing in details the 
mass matrices obtained in such a situation, it is useful to consider the relevant corrections 
coming from three-flavon insertion operators and from the vev shifts quoted in eq. 02Up . The 
relevant feature of such corrections is the filling of the vanishing element of the mass matrix 
in eq. flMl) . which turns out to be se^, that is of Oi^X^). This allows to shift 6*23 away from 
zero by a quantity of O(A^). The relevant operators are of the form FT^H^ {^i^l + Xi^l) 
computed with leading order vevs and FT^H^ (v^^v^^ + '{'vXi) with flavon vevs at the next to 
leading order (see appendix B for details). Putting together all these elements we are now 
in the position to write down the most general mass matrix allowed in our model: 

A^ A^ A^ \ 
rrid^VdX'' I A3 A2 A2 . (29) 

A A2 1 / 

The unitary left-handed rotation Ud is obtained diagonalizing nidrrid whereas the right- 
handed one Ui is the charged lepton rotation. In terms of the Cabibbo angle they are given 
by: 



[/d ~ A2 1 AM f/^ ~ A 1 A2 . (30) 





We see that Ud perfectly reproduces the correct order of magnitude of Vub and Vcb, whereas 
Vus turns out to be a bit smaller than expected, as in many other models based on non- 
abelian discrete symmetries. On the other hand, in the construction of the neutrino mixing 
matrix Upmns [20], Ui will induce corrections of 0{X) to 9^ and 6*13 and O(A^) to 623. These 
corrections turns out to be relevant to shift the solar and reactor angles from their BM values 
to the experimental ones. 

For the sake of completeness, we recompute the previous quantities including all C(l) 
coefficients. The mass matrix of the down quarks is now given by: 



rrid ~ ctb Vd A^ 



Xi\^ X2A^ x^X^ 

X/^\^ xsA^ —x^\^ 

Xq\ XjX^ 1 



(31) 



where ab is the bottom Yukawa coupling from the first operator in eq. ( 12T]) and Xi are linear 
combinations of the other Yukawa couplings of eq.( l2T]) . rescaled by ab- 

Xl = -71 + 72 - 73 - 75 - 76 - 7ll + 7l2 - 7l3 + IS^ui + 74^9=1 
X2 = -78 
X3 = 74 + 78 

X, = -/34 + /35 - /36 + /3Ai (32) 

X5 = -/3i 

xq = — «! — ^2 + ab5^^ 

2x7 = 3a5 + VSar - 2ai{6uj^ + 6^J + 2a2{Suj^ + 6^J . 

Notice that the (22) and (23) elements are opposite to each other; this equality is broken 
only at C(A^), for example, by operators of the form FT26 {ipj)2fe- From eq.( l3T]) it is easy 
to obtain the expression for the charged lepton mass matrix by transposition and changing 
accordingly xi — X5 with x'l — x'^, which differ from the previous ones by 0{1) coefficients. 
This is the effect of introducing the copies T[ 2 of the first two tenplet fields, whose zero modes 
are different from those of Ti 2 and couple with the charged leptons only. Summarizing, the 
mass matrix of the charged leptons is: 



rrie ~ ab Vd A^ 
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(33) 



Working for simplicity in the limit of real coefficients, the quark and charged lepton masses 
in unit of a^ Vd A^ are explicitly given by: 
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We see that, at the GUT scale, the h — r unification is recovered (up to very small corrections 
of (9(A*^) not listed here); also, the other two mass ratios md/rrie and ms/m^ are both of the 
same order of magnitude but not strictly equal because of different 0{1) coefficients. 



Keeping only the leading order for each matrix elements, the left-handed rotation Ui for 
charged leptons is given by: 
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4 Up quarks mass matrix 



The up quark mass matrix is completely determined by operators with no more than two- 
fiavon insertions. In the following we list the relevant non-vanishing terms in the superpo- 
tential disregarding those fiavon contractions which will give contributions only at the NLO, 
not relevant for our discussion: 
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Taking into account the conditions in eq. (l28|l . the corresponding mass matrix reads: 



m. 
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(r/i + r72)A8 (ri + r2)A6 aX^ 
vlatX\ (ri + r2)A6 pX^ -^A^ 

aX^ ~6X^ 1 



(38) 



where, once again, we have rescaled the Lagrangian coefficients by at and used the same 
symbols. It is easy now to read masses and mixing matrix; for the first we have: 



nit : rric : ^u ~ 1 : A"^ : A^ 



(39) 



which perfectly matches with the GUT expectations. To avoid large dimensionless coeffi- 
cients, we assume here that Vu,d ~ Af ° ^. We are allowed to do that because of the freedom 
related to the boundary values f°^; in fact, the electroweak scale is determined by the 
relations: 
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and, unless {H^K,{y)) are constant in the fifth coordinate, one would expect wj]^ 7^ Vu,d- 
With this assumption, the Yukawa coupling of the top quark is of order one and, thanks 



to eqs. 

diagonalizing mupml^ is 



, also all the other couplings are of the same order. Up to 0{\ ) the matrix 
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(41) 



We see that the up sector contributes to both {Vus — Kd) and {Vub — Vtd) of the CKM, which 
in turn results to be: 
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It is interesting to observe that the different coefficients of the (13) and (31) elements 
can explain the experimental difference among these matrix elemental. As a final comment 
of this section, we observe that: 
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(43) 



which can easily reproduce the experimental value mb/rrit ~ A^ for v^/v^ ~ 1/A (and 0{1) 
Yukawa couplings). Considering the scaling Vu,d ~ -^"^ud) the condition v^/v^ ~ 1/A implies 
tan/3 ~ 5. 

5 The neutrino sector 

In this section we show how to get a description of the neutrino masses and mixings based 
on the 5*4 symmetry. As already outlined in the Introduction, the mixing matrix at leading 
order will have the BM structure. Our main issue here is to show that dimension 5 Weinberg 
operators are enough to get the BM mixing matrix at LO and a mass spectrum compatible 
with the data. Since a see-saw version (with subdominant D=5 contributions) has been 
carefully studied in pTl], in the second part of this section we limit ourselves to show that 
such a mechanism can also be successfully employed in the GUT version of the model, with 
identical LO results. 

5.1 The neutrino sector from effective operators only 

We start writing the leading and next to leading order contributions to neutrino masses from 
the Weinberg operator: 



w^ = ^^{FF),H,H, + y^{FF),,H,H,ip, + y^(FF),H,H,^, 



(44) 



^In addition, to reproduce the correct Cabibbo angle with 0(1) coefficients one has to ask, for example, 
for a negative X2 or x^ . In this case, it is enough to take the absolute value of these coefficients equal to the 
unity to gain an enhancement of a factor of 3. 
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from which the mass matrix is (we use the symbol A according to eqs. 



): 



m„ 



A 



z/«,iA yw + yw2^ 

y«,iA yw + yw2^ 



(45) 



This matrix is exactly diagonalized by BM and the eigenvalues are (replacing sv^ with v^) 



nil 



m2 



1713 



+ iyw2 - V2y^^)\ 
+ {yw2 + V2y^^)\ 



(46) 



-f ] {yw + yw2^)- 



Notice that they satisfy the sum-rule: 

mi + 1712 + 2171-3 = , 



(47) 



which allows to reduce the number of independent parameters in the neutrino mass matrix. 
We can reparametrize the masses in terms of two different complex parameters: 



V2y^^X = B = be'^'' 



a,b> 



(4J 



from which: 



\mi\ 



\m2\ 



ms 



(2 \ 
^j {a^ + b^ - 2 ab COS A) 

,2X2 

^\ {a^ + b^ + 2 ab cos A) 



A 

„,2\ 2 



(49) 



A' "' 



where A 



%. Since the solar mass difference is given by: 

,2\ 2 



Am: 



ol 



\m2\ 



niiY = 4: [-^\ ab cos A 



(50) 



the condition cos A > must be fulfilled. This implies that the spectrum is of inverted type 
because the condition Im^l > \m2\ cannot be satisfied. We then define the quantity 



^^Lm 



,2\ 2 



m^l^-lmar = ( -^ ) 6(6-2acosA) 



(51) 



which is positive for 



cos A < 



2a 



The value of r 



4a cos A 



AmL 

Am\^.^ b — 2a cos A 



(52) 



(53) 
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can be made small only if cos A ~ because the other possible conditions b » 2a cos A 
reduces to small cos A and h « 2a cos A is in conflict with eq. fl52|) . We do not consider 
here the possibility of a << 1 because this parameter is dominated by yyj and we prefer to 
work with 0{1) coefficients in the superpotential. 
The experimental value r ~ 1/30 is reproduced if 



cos A 



2a(2 + r; 



r ~ 10 



-2 



(54) 



for a ~ h ~ 0{1). Notice that, using this relation into eq.f lSTj) . we can estimate the order of 
magnitude of the scale A, which turns out to be: 



A 



hvl 



a/ Am, 



2 
atm 



lO^^GeV^ 



(55) 



for Vu = 100 GeV and Am^^„ = 2.4 x 10"^ eV^. 



At the next to leading order, we have to consider two-flavon insertion operators and the 
shift of the flavon (p^, of eq. flT6|) . that is 6^^^. The contributions of the new operators to the 
superpotential is: 



6w,, 



^"^ {FF),{^,^,),H,H, + ^{FFU^,^,)2H,H, + ^{FF),,{^,^,h,H,H, + 



A4 

Uwd 



A4 



A4 



(56) 



Vwj 



— {FFUtH.H, + !^{FF^M.H,H, , 



with the corresponding mass matrix given by 






7A 






-y 



«)i 



A 



y 



Wj 



X' 



yy, + yy,^\ + y^y 



(57) 



1(32/^4 + 2?/^j5^JA2 

where y'^ = -2y^^ + y^^ + y^^ and yl = l/2{-Ay.^^ - y^^ + 2y^^). It is easy to show that, 
in the limit 6^^ —^0, the mass matrix is still diagonalized by BM. For 6^^ 7^ 0, we observe 
that 61^ and 612 are not corrected whereas ^23 is shifted by a quantity of (9(A^): 



^23 



n 

4 + V2y, 



ywi 



6uiX 



(5^ 



To get the neutrino mixing matrix, we also have to take into account the corrections coming 
from the charged leptons, eq. fl35|) : to be consistent up to C(A^), we have to include all C(A^) 
terms in eq. (l35l) . which mainly affect the diagonal elements (the elements of 0{X) are not 
modified at this order). Since such corrections are important for 623, what matter for us are 
the {1/1)22 and {Ui)s3 elements, which turn out to be {U()22 = 1 — (2^4 + ^'^^q)"^ / i2x'^) A^ and 
(^^)33 = 1 ~ (^^e/^) ^^- With this in mind, the neutrino mixing angles can be written as: 
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13 



723 



712 



1 

72 



X 



X 



A 



vr 



+ 






ywx 

2yw 



(59) 



IT 

4 



1 

V2 



x'^ + 2 



4^6 \ ^ 
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We explicitly realized the pattern in eq.(|2]), with all corrections computed analytically. These 
relations show a certain amount of correlation; for example, the solar angle can be written 
as: 



9 



12 



IT 

4 



^13 - \/2xqX 



(60) 



so that we expect a large deviation of 612 from maximal mixing when ^13 is significantly 
different from zeros, unless some cancellations with the xq parameter is at work. This is 
confirmed by a numerical analysis by treating the parameters in eqs.f l35p and ( 157|) as random 
complex numbers. Since it is crucial that the corrections from the charged lepton sectors 
have coefficients of 0{1), we restrict the absolute values of the Xj between 1/2 and 2; all the 
other parameters are extracted in the range of absolute values between and 3. Notice that 
the results in eqs. (J35l) and ( l57l) are valid at the GUT scale; however, since the model predicts 
inverted hierarchy, the effect of the running to the SUSY scale can be safely neglected. The 
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Figure 1: sin^6'i2 as a function of sin^ 9 13 (left panel) and siia^ 623 vs sin^6'i3 (right panel). 
The gray bands represents the regions excluded by the experimental data at 2a on sin^ Q\2 
and sin^ ^23 \^ whereas the vertical lines enclose the 2a range on siv? 6*13 as obtained in J^. 
See text for further details. 

results of such an analysis are shown in Fig.([T]); in the left panel, we plot sin^6'i2 as a 
function of sin^ ^13. The gray bands represent the regions excluded by the experimental data 
at 2(j on sin^ 612 [S] whereas the vertical lines enclose the 2a range on sin^ ^13 as obtained 
in [5] (old reactor fluxes). The points correspond to choices of the parameters reproducing 
Am^j^, Am^g; and the mixing angles within a 3cr interval, for the inverted ordering of the 
neutrino mass spectrum. We see that a relevant fraction of points fall in the box allowed by 
the current limits on neutrino parameters. In the right plot, we show sin^ 6'23 (at 2a from 
[5]) as a function of sin^6'i3. To guide the eye, we also display an interval of ±A^ around 
sin^ ^23 = 1/2. These results confirm the analytical formulae in eq. (l59|) : within this model 
it is difficult to produce large deviations of ^23 from maximal mixing since large corrections 
can only be produced invoking some fine-tuning among 5^^ and the ?/j's. 

In Fig.(j2]) we plot the model prediction for |mee| as a function of the absolute value 
of the lightest mass m^. The shaded area corresponds to the region allowed by current 
neutrino data, for a mass ordering of inverted type. The vertical band corresponds to the 
future sensitivity on the lightest neutrino mass of 0.2 eV from the KATRIN experiment 
[2T] and the horizontal line to the future sensitivity of 15 meV of CUORE [22]. The figure 
suggests that many of the model prediction points are in the quasi-degenerate region, above 
the sensitivity of CUORE. Furthermore, the values of the lightest neutrino mass rely in a 



13 



region close to 



2 
atm 



mal ~ A/Aml^, as it can be understood from eq. fH^ and recalling that 
As a consequence, much smaller values of |?7i3| can only be obtained if, 
thanks to additional fine tunings of the parameters, a cancellation between the NLO and LO 
contributions takes place. Having obtained this only for very few points, we can consider 
the indication for a lower bound on \m^\ around 5 ■ 10~^ eV as a reasonable prediction of the 
model. 
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Figure 2: |?7iee| os o function of the lightest mass m^. The shaded area corresponds to the 
region allowed by current neutrino data, for a mass ordering of inverted type. The vertical 
band corresponds to the future sensitivity on the lightest neutrino mass of 0.2 eV from the 
KATRIN experiment WT^ and the horizontal line to the future sensitivity of 15 meV of 
CUORE 



5.2 The see-saw version 

We introduce three right-handed neutrinos transforming under SU(5)xS'4 x Z^ as (l,3i, 1) 
and carrying U(1)r = +1. At leading order the superpotential reads as follows: 



w,. 



AV2 



Fu'H^ + M.Ai/V^ + a{u'u%,ifi, + biu'u'^U, 



(61) 



and the resulting Dirac and Majorana mass matrices are: 



m 



D 



syvl 




M: 



N 



M + bX -aX -aX 

-aX M + bX 

-aX M + bX 



A 



(62) 



similar to those quoted in [TT]. In particular, the expression of the eigenvalues are the same 
and the spectrum results compatible with normal hierarchy only. We refer to that paper for 
a detailed discussion of the leading order results implied by eq.(l62l) (and, in particular, of a 
proof that D=5 operators can be made subdominant). At the next level of approximation. 
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the additional contributions to w^, are: 



Sw,, 






(63) 



+f((^=^^ 



ii^^^u)^) + -^{.^''^''Vu)iu + h.c. 



Since there is no mixing between the flavon fields of different sectors, these contributions 
only modify the expressions of the neutrino masses but not the mixing matrix, which is still 
of BM form. However, as shown in eq. (l2UI) . the fiavon Lp^ takes a non- vanishing component 
in the first element, which is the source of the BM breaking. This affects the third operator 
in eq.f l6T|) and the first operator in eq.f l63p and it is of the same order of magnitude in A of 
the contributions of the operators in eq.f l63|) with coefficients ai, ..,04. To be more explicit 
and put in evidence the sources of the BM breaking, we write the leading deviation of the 
light neutrino mass matrix in the following way: 



m. 



w 



X 


y y 


y 


Z ~\~ -/I X Z 


y 


X — z z — A 



(64) 



where 



A 



yy 

M2A 



{2M,yi - ay,) 6,, X' 



(65) 



and X, y, z are linear combinations of the superpotential parameters. It is clear that the 
limit y4 — 7- reproduces a mass matrix diagonalized by BM. A straightforward calculation 
of the eigenvectors up to (!?(A^) shows that the solar and reactor angles do not receive any 
corrections whereas 623 deviates from maximal mixing by a quantity proportional to 6,-^^ A^. 
Taking into account the corrections from the charged leptons, we arrive at a final result 
similar to eq. (l59l) : 



e 



13 



723 



'12 



1 
71 



TT 

4 + 



Xa 



Jb c; 



A 



X A 



4x'2 2 



Xj 



'y\ 



2M„ 



yi 

yv 



X' 



(66) 



TT 

4 



1 

V2 



X4 + 2 X5X6 \ 

4 r- 



6 Conclusions 

In this paper we have constructed an SU(5) model based on the fiavour symmetry S4X Z^x 
U{1)fn, where the BM mixing is realized at the LO in a natural way. In order to get a 
realistic pattern of fermion masses and mixing and to get rid of the usual problems in SU(5) 
constructions, like the rigid mass relation rrie = rn^ and the doublet-triplet splitting problem, 
we have embedded the model in five dimensions, with the fifth dimension compactified on 
the orbifold S^/Z2. In this way, we were able to reproduce the correct mass ratios in both 
charged lepton and quark sectors as well as a good agreement with the entries of the CKM 
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mixing matrix, with only a moderate fine-tuning needed to describe the Cabibbo angle. Since 
exact BM mixing is not confirmed by neutrino oscillation data, large corrections are needed 
in order to lower the value of the solar angle. Such corrections arise at NLO through the 
diagonalization of the charged lepton mass matrix while, at the same time, the reactor angle 
receives corrections of 0{Xc), which are perfectly in agreement with a large ^13 emerging 
from global fits to the world neutrino data. An important feature of our model is that the 
shift of sin^6'23 from the maximal mixing value of 1/2 is expected to be of 0{Xq) at most. 
In order to reproduce the experimental value of the small parameter r = Am'j.^i/Am'^^,^ we 
need some amount of fine tuning. Here we choose to tune the phases of two complex numbers 
entering the neutrino mass matrix at leading order; this implies that the neutrino spectrum 
is mainly of the inverted hierarchy type (or moderately degenerate) and that the smallest 
light neutrino mass and the Oz//3/3 parameter |mee| are expected to be larger than about 0.01 
eV, as we obtained from a numerical analysis with mass matrices at the NLO. 
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Appendix A: corrections to the leading order vacuum 
alignment 

In this appendix we explicitly compute the correction to the vacuum alignment of eq.(l9]). 
For the sake of completeness, we use the following parametrization of the fiavon vevs: 
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The superpotential responsible for the vev shifts reads as follows: 

a7ix'i)2Xi + asixDsiXi + a9^Avl)3i + aioi^l)2Xi + «ii(v5^)3iX^+ 

C [(^li^Ds^Vu + P2il + HvlKv, + Mu{vl)i + PMhM + (67) 

^° [ll{^])-ii^v + l2iu{^])2 + lAxVl^^^u + lAiu{xl^)2 + 
75(V5vV5£)3iX^ + 76(V5i.V5£)32X£ + lliuiV 1X1)2] + 



hi'^u'^dlXl + '^6(V'!.V'£)2X£ + ^li'^u'^d^iXt + '^8(V^i/V5^)32X^ + 59^i/(V^^X^) 



32 J 



Even after considering the leading order relations in eq. fllOl) . some of the minimizing equa- 
tions are still cumbersome. In the following we concentrate on those which easily allow to 
extract the relevant information on the vev shifts: 



dAw, 



d 



2^i<(^.2+'^.3)=0 






(68) 


^x 


^f?,Vvi{K3+^X2) - 


- 2/25x2^X 


- Vi^^A, 


(69) 


^v^ 


,v^{5^,-5^^) = Q. 






(70) 



dAwd 
dAwd 

We see that the last equation implies ^^3 = 5^^^ which, using eq.(l69l). forces 5^^ = S^^ = 0- 
The first equation tells us that 6,^2 = ~^i^3 ^^^ ^^at they are undetermined. Using the 
previous relations, the remaining minimizing equations allow to fix the magnitude of the 
other vev shifts in terms of the superpotential parameters. In particular, from dAwd/dip^^ 
and dAwd/d(f^,^ we get a set of conditions on 6^-^ and 6v^, which result both different from 
zero. From dAwd/dx2 3 = we obtain an expression for 6^-^ and 6^-^ , respectively; the 
last equation dAwd/ipe^ = depends on 5^2, which is also completely determined by the 
superpotential parameters. 



17 



Appendix B: corrections to the fermion mass matrices 

In this appendix we quote all the operators with three-flavon insertions responsible for the 
down quark sector: 

FTsH^ : (<^^)i<^^, {^l)2^e, {y^Ds^^e, {^l)2Xe 

FT2H^e : (<^^)i<^v, {^l)2^v, {^Ds.ifiu, {^^i)i^e 

{f'i)2fi, (¥'^)3i^£, {xl)2Xi, ixDsiXh {'pDsi^u, Cfu, 

{v^e)2Xi, {vl)3iX£, {xDiVe., {x^e)2Vi, {x^e)3i¥>i 

FT2Hid' : {vi)iVu, {vi)2Vy. ivDs.V,., ivDsi^iy, {xDiVi^, {xl)2<^u, 

ixDsi^u, ix^i)3i^u, {^uVi)2Xe, {VuVe)3iXi, 

(71) 

{^i^,y)32Xe, {vae)3^iu 

FTiH^e^ : the same as FT2H-^e' 

FTiH-^e'^e' : the same as FT3H-5 

FTiH-^99''^ : the same as FTgi/g^ 

FTiH^e'^ : the same as FTai/g^' 

Appendix C: the group ^4 

We report here the multiplication table for 5*4 and we list the Clebsch-Gordan coefficients 
in the basis used in the paper [11]. In the following we use a^ to indicate the elements of 
the first representation of the product and (3i to indicate those of the second representation. 
Explicit forms of S and T are as follows. In the representation 1 we have T = 1 and 5* = 1, 
while T = — 1 and S* = — 1 in 1'. In the representation 2 we have: 



For the representation 3, the generators are: 

[ ' -71 -A \ 

(73) 



-1 0\ / -71 -72^ 



T= I -i I S 

i ^ 
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i_ 1 _i 

v^ 2 2 

\ V2 2 2 ) 



In the representation 3' the generators S and T are simply opposite in sign with respect to 
those in the 3. 

We start with all the multiplication rules which include the 1-dimensional representations: 

1 Rep = Rep ^ 1 = Rep with Rep any representation 



h 
h 



af3 



a (5 2 
—af3i 



>l2 = l 
) 3i = 3^ 



I2 ® 32 = 3i ^ 

The multiplication rules with the 2-dimensional representation are the following ones: 

' 1 ~ ai/3i + 0:2/32 
aij32 - a2l3i 




l©lc 



with < 



2®3i = 3i 



2 ® 32 = 3i © 32 



with 



with < 



I2 

2- 



02/32 — cei(3i 
ai(32 + a2(3i 



Ol/?! 

^02/33 - |tti/32 



^02/32 - ittl/33 



-a2/3i 
^ai/33 + |a2/32 
^ai/32 + 1^2/33 

-a2(3i 
^ai/^s + |a2/32 
^ai/32 + |«2/33 



«i/3i 
1 

2' 



V3 
2 

^02/32 - haiPs 
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The multiplication rules involving the 3-dimensional representations are: 

1 ~ ai/3i + 0:2/33 + "3/92 

ttlA - K«2/33 + «3/32) 
^(02/52 + a3/33) 



3i®3i 



l©2©3i 



with < 



"3/33 - "2/32 

3i ~ ( aiPs + as/?! 
-tti/32 - tt2/3l 

a3/32 - a2/33 
tt2/3i - ttl/32 
tti/33 - "3/3l 



3i®32 = l2©2©3i©32 



L2 ~ ai(3i + 02/33 + ^3/32 



V3 



with < 



_ («2/32 + ^3/33) 
-ttl/3l + |(«2/33 + «3/32 



03/32 - O2/33 

3i ~ ( 02/3i - oi/32 

ttl/33 - "3/31 
"3/33 - "2/32 

"1/33 + "3/3l 

— Oi/32 — 02/3l 
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